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Parametrization of nuclear parton distributions is investigated in the leading order of a s . 
The parton distributions are provided at Q 2 —l GeV 2 with a number of parameters, which are 
determined by a x 2 analysis of the data on nuclear structure functions. Quadratic or cubic 
functional form is assumed for the initial distributions. Although valence quark distributions in 
the medium x region are relatively well determined, the small x distributions depend slightly on 
the assumed functional form. It is difficult to determine the antiquark distributions at medium 
x and gluon distributions. From the analysis, we propose parton distributions at Q 2 —l GeV 2 
for nuclei from deuteron to heavy ones with the mass number A ~ 208. They are provided 
either analytical expressions or computer subroutines for practical usage. Our studies should 
be important for understanding the physics mechanism of the nuclear modification and also 
for applications to heavy-ion reactions. This kind of nuclear parametrization should also affect 
existing parametrization studies in the nucleon because "nuclear" data are partially used for 
obtaining the optimum distributions in the "nucleon" . 
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I. INTRODUCTION 

Unpolarized parton distributions in the nucleon are 
now well determined in the region from very small x 
to large x by using various experimental data. There 
are abundant data on electron and muon deep inelastic 
scattering. In addition, there are available data from 
neutrino reactions, Drell-Yan processes, W productions, 
direct photon productions, and others. The optimum 
parton distributions are determined so as to fit these ex- 
perimental data. Initial distributions are assumed at a 
fixed Q 2 with parameters, which are determined by a 
X 2 analysis. Now, there are three major groups, CTEQ 
(Coordinated Theoretical/Experimental Project on QCD 
Phenomenology and Tests of the Standard Model) [1], 
GRV (Gliick, Reya, and Vogt) [2], and MRS (Martin, 
Roberts, and Stirling) [3], which have been investigating 
the unpolarized parametrization. 

Because the distributions themselves are associ- 
ated with nonperturbative Quantum Chromodynamics 
(QCD), they cannot be described precisely by theoretical 
methods at this stage. Therefore, the determination en- 
ables us to understand internal structure of the nucleon, 
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hence to test the nonperturbative hadron models. In ad- 
dition, the parton distributions are always necessary for 
calculating cross sections of high-energy reactions involv- 
ing a nucleon. Furthermore, there is other importance 
in determining future physics direction. If the distri- 
butions are precisely known, any experimental deviation 
from theoretical predictions should indicate new physics 
beyond the standard model, for example a signature of 
subquark system. 

The situation is worse in polarized distributions. Lon- 
gitudinally polarized distributions have been investigated 
in the last decade. From the analyses of many polar- 
ized electron and muon deep inelastic experimental data, 
several parametrizations have been proposed. Neutrino, 
Drell-Yan, and other data are not available unlike the 
unpolarized studies, so that antiquark and gluon distri- 
butions cannot be determined accurately at this stage [4] . 
More precise determination should be done by the polar- 
ized Relativistic Heavy Ion Collider (RHIC) experiments 
in the near future. 

It is well known that nuclear parton distributions are 
modified from the corresponding ones in the nucleon 
according to the measurements of nuclear F2 structure 
functions [5] . Because F2 is expected to be dominated by 
sea- and valence-quark distributions at small and large x, 
respectively, the modification of each distribution should 
be determined by the F2 measurements in the small or 
large x region. However, it is not straightforward to find 
the detailed x dependence of these distributions. Fur- 
thermore, it is also not obvious how nuclear gluon distri- 
butions are constrained by the F2 data. There are tri- 
als to obtain the distributions from the data in a model 
dependent way [6] and in a model independent way by 
Eskola, Kolhinen, and Ruuskanen [7]; however, they are 
not x 2 analyses. Therefore, this paper is intended to pio- 
neer the x 2 analysis study for obtaining optimum nuclear 
parton distributions [8]. 

There are following motivations for investigating the 
nuclear parametrization. The first purpose is to test var- 
ious nuclear models for describing the nuclear structure 
functions. From the detailed comparison, the most ap- 
propriate nuclear model could be determined in the high- 
energy region. Furthermore, we may be able to find an 
explicit subnucleon signature in nuclear physics. The sec- 
ond purpose is to calculate cross sections of high-energy 
nuclear reactions accurately. For example, precise initial 
conditions must be known in heavy-ion reactions for find- 
ing a signature of quark-gluon plasma. In many calcu- 
lations, nuclear parton distributions are simply assumed 
to be the same as the corresponding ones in the nucleon. 
The third purpose is related to the aforementioned nu- 
cleon parametrization. In obtaining parton distributions 
in the "nucleon" , some "nuclear" data are actually used 
without considering the nuclear modification. For ex- 
ample, neutrino F3 data are essential for determining 
the valence-quark distributions in the nucleon. However, 
most data are taken for the iron target! In future, ac- 
curate neutrino-proton scattering data could be taken if 



a neutrino factory is realized [9]. It is inevitable to uti- 
lize the nuclear data at this stage. Therefore, our nu- 
clear studies should be useful for improving the present 
parametrizations in the nucleon. 

This paper consists of the following. First, dependence 
of the mass number A and scaling variable x is discussed 
in II. Then, our analysis method is explained in III, and 
results are presented in IV. Obtained distributions are 
provided in Sec. V for practical usage. Finally, our nu- 
clear parametrization studies are summarized in section 
VI. 



II. FUNCTIONAL FORM OF NUCLEAR 
PARTON DISTRIBUTIONS 

Because there is no prior x 2 analysis on the nuclear 
parametrization, we should inevitably take the A and x 
dependence as simple as possible for the first step trial. 
However, the functional form of x should be taken inde- 
pendently from any theoretical nuclear models as a fair 
analysis. Our standpoint is to test the models without 
relying on them. We discuss an appropriate functional 
form of the distributions in the following subsections. 



A. A dependence 

There is some consensus on physics mechanism of nu- 
clear modification in each x region. Because different 
mechanisms contribute to the modification depending on 
the x region, the A dependence could vary in different x 
regions. However, in order to simplify the analysis, we 
introduce a rather bold assumption: The A dependence 
is assumed to be proportional to 1/A 1 ^ 3 . Physics behind 
this idea is discussed by Day and Sick in Ref. [10]. In any 
nuclear reaction, the cross section is expressed in terms 
of nuclear volume and surface contributions: 

a A = Aa v +A 2 / 3 a s . (2.1) 

Therefore, the cross section per nucleon is given as 



A 



: cry + 



1 



A 1 / 3 



(2.2) 



If cry and as depend weakly on A, the 1/A 1 / 3 depen- 
dence makes sense as the leading approximation. In fact, 
according to the measured F^ / F 2 D data, this A depen- 
dence is justified [10]. 

In order to illustrate the justification, we show actual 
F 2/ F 2 data as a function of 1 - 1/A 1/3 at x = 0.5 to- 
gether with a straight line in Fig. 1. It is obvious from 
the figure that measured ratios are well reproduced by the 
line with the 1 — 1/A 1 / 3 dependence. We also looked at 
other regions, for example, the small x region and found 
that the data could be described by this A dependence. 
In this way, we decided to employ this simple assumption. 
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FIG. 1: A dependence of measured F<^/F<f ratios at x — 0.5. 



Of course, the 1/A 1 / 3 dependence may be an oversimpli- 
fication. Detailed studies of the A dependence are left as 
a future research topic. 

We should mention, however, that the exact 1 — 1/A 1 / 3 
dependence is not completely consistent with the condi- 
tions of nuclear charge, baryon number, and momentum. 
We need to assign fine-tuning parameters which adjust 
the nuclear dependence. The details will be explained in 
the next subsection below Eq. (2.8). 



B. x dependence 

Available nuclear data on the structure function F£ 
are taken in fixed-target experiments at this stage, and 
they are shown as a function of Q 2 and x = Q 2 /(2raz/), 
where v is the energy transfer, m is the nucleon mass, 
and Q 2 is given by Q 2 = — q 2 with the virtual photon 
momentum q. The initial nuclear parton distributions 
are provided at a fixed Q 2 (= Qq), and they are taken as 



f l A (x,Q 2 ) = w l (x,A,Z)f i (x,Ql), 



(2.3) 



where Qq) is the parton distribution with type i in 

the nucleus A and /i(^,Qo) * s the corresponding parton 
distribution in the nucleon. We call Wi(x, A, Z) a weight 
function, which takes into account the nuclear modifica- 
tion. 

We lose a piece of information by using Eq. (2.3). A 
scaling variable for the lepton-nucleus scattering could be 
given by xa = Q 2 /(2pA • q), where pa is the momentum 
of the target nucleus. We could define another variable 
for lepton- nucleon scattering as xjy = Q 2 /(2p • q) with 
the nucleon momentum p within the nucleus. The vari- 
able xa indicates the momentum fraction of a struck- 
quark in the nucleus and it is kinematically restricted 
as < xa < 1. Because of these definitions, we have 
the relation x = xaMa/tu, where Ma is the nuclear 
mass, in a fixed target reaction. In this way, we find 
that the range of nuclear parton distributions is given 
by < x < A. Therefore, the extremely large x region 



(x > 1) cannot be described in our present approach with 
Eq. (2.3) because the distributions in the nucleon vanish: 
/i(x>l,Q§) = 0. 

There are following reasons for using Eq. (2.3) irre- 
spective of this issue. First, there are not so many exper- 
imental data in this x region. Second, even if finite parton 
distributions are assumed at the initial point Qq, there 
is no reliable theoretical tool to evolve them to larger Q 2 
points where the data exist. Therefore, such large x dis- 
tributions cannot be accommodated in our x 2 analysis. 
There is a nuclear Q 2 evolution code in Ref. [11]; how- 
ever, the original evolution equations are written in the 
range < x < 1 [12]. Third, because the nuclear modifi- 
cation is generally in the 10—30% range for medium and 
large size nuclei, it is much easier to investigate the modi- 
fication from the distributions in the nucleon rather than 
the absolute nuclear distributions themselves. Forth, the 
structure functions F^ or parton distributions are tiny 
at x > 1, so that they do not affect the calculated cross 
sections significantly unless an extreme kinematical con- 
dition is chosen. 

Now, we proceed to the actual x dependence. 

(1) Nuclear modification oc 1 — 1/A 1 / 3 

As concluded in the previous subsection, the nuclear 
modification part of Wi(x, A, Z) is assumed to be propor- 
tional to 1 - 1/A 1 / 3 : 

Wi(x, A, Z) = 1 + ^1 — (function of x, A, and Z). 

(2.4) 

(2) Introduction of 1/(1 — x)^ factor 

The nuclear parton distributions have finite values 
even at x = 1 in principle, so that the weight function 
should behave as 



Wi(x 1,A,Z) 



(2.5) 



This equation should hold whatever the modification 
mechanism is. In order to reproduce this feature, the 
factor 1/(1 — x)P* is introduced in the x dependent func- 
tion part of Eq. (2.4) with a parameter /%(> 0): 



(function of x, A, and Z) oc 



(l-x) 



ft ' 



(2.6) 



(3) Saturation of shadowing or antishadowing 

We assume saturation of the function Wi(x, A, Z) at 
x — > 0. It is considered to be a reasonable assumption 
unless there is a peculiar mechanism to produce singular 
behavior at small x. As far as the experimental F^/F® 
data suggest, the shadowing at small x tends to satu- 
rate as x — > 0. If this saturation is assumed, the weight 
function w* becomes 



Wi(x^0,A,Z) = l+[l- 



A 1 / 3 



ai(A,Z), (2.7) 



where ai is the parameter which controls shadowing or 
antishadowing magnitude. 
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From these discussions, an appropriate functional form 
becomes 



Wi(x,A, Z) = 1+1 



1 \ Oi(A, Z) + Hi(x) 



A 1 / 3 



(1 - x)P< 



(2.8) 



where the additional x dependent part is written as 
Hi(x), and it determines minute x dependent shape. The 
reason why we let the parameters ai depend on A and 
Z is the following. As explained in the next section, 
there are three obvious constraints, charge, baryon num- 
ber, and total momentum, on the nuclear distributions. 
If these conditions are satisfied for a nucleus, they are 
also fulfilled for other nuclei with the same Zj A ratio. 
For example, the conditions could be satisfied for all the 
isoscalar nuclei. However, we also analyze nuclei with 
different Zj A ratios. Then, even if the three conditions 
are satisfied in a certain nucleus, they are not strictly 
fulfilled for other nuclei with different Z/A factors. In 
order to avoid this kind of failure, nuclear dependence is 
assumed for the parameters ai in the valence-quark and 
gluon functions. They are determined by the x 2 analysis 
so as to satisfy the three conditions for any nucleus. 

Because this is the first x 2 trial, we would like to 
simplify the functional form of Hi(x) as much as pos- 
sible. A simple idea is to expand it in a polynomial form 

Hi(x) = biX + CiX 2 H . An advantage of this functional 

form is that the polynomials of x and 1 — x are very 
easy to be handled in the Mellin-transformation method 
of the Q 2 evolution. Because a direct x-space solution 
[11] is used in our Q 2 evolution, it does not matter in 
our present analysis. However, it is important for public 
usage if we consider the fact that many researchers use 
the Mellin transformation as their Q 2 evolution method. 

It is obvious that Hi(x) = biX cannot explain the com- 
plicated F^/F® shape in the medium x region. There- 
fore, the simplest yet realistic choice is to take Hi(x) = 
biX + CiX 2 . This seems to be acceptable in the sense that 
the medium x depletion of F^ / F® can be explained to- 
gether with the F<2 shadowing at small x. However, if the 
depletion is described by the valence-quark behavior as 
our common sense suggests, the valence-quark distribu- 
tions show anti-shadowing at small x due to the baryon- 
number conservation. As explained in Refs.[9, 13], it is 
not obvious at this stage whether the valence distribu- 
tions indicate shadowing or antishadowing without accu- 
rate neutrino-deuteron scattering data and also small x 
nuclear data. Therefore, this functional form should be 
considered as one of possible options: 



Wi(x,A,Z) = 1+1 



A 1 / 3 
"quadratic type" . 



a,i(A, Z) + biX + CiX 2 



(1 - x)P* 



(2.9) 



We call this weight function the quadratic type in the 
following discussions. 

If the next polynomial term diX 3 is added in addition, 



the function becomes more flexible in fitting the data: 

1 \ di(A, Z) + biX + ^x 2 + diX 3 



Wi(x,A,Z) = 1+1 



A 1 / 3 



(1 - x)& 



"cubic type" . 



(2.10) 



We call this weight function the cubic type. An ad- 
vantage of the additional term is, for example, that the 
weight function becomes flexible enough to accommodate 
both possibilities, shadowing and anti-shadowing, in the 
valence-quark distributions. A disadvantage is that the 
number of total parameters becomes larger, so that the 
total computing time becomes longer. 

These quadratic and cubic functional types are used 
in our x 2 analyses. Finding the optimum point for the 
parametrization set, we expect to explain the major fea- 
tures of the measured F 2 A / F® ratios. 



III. ANALYSIS METHOD 

In addition to the initial functional form, there are 
other important factors for performing the y 2 analysis. 
In this section, the details are discussed on used experi- 
mental data and our x 2 analysis method. 



A. Experimental data 

There are many available experimental data which 
could indicate nuclear parton distributions. However, 
we restrict the data to those taken by the deep inelastic 
electron and muon scattering. Neutrino, Drell-Yan, and 
other hadron-collider data are not used in our present 
analysis with the following reasons. At first, we would 
like to investigate how the distributions are determined 
solely by the electron and muon scattering data. Next, as 
it is mentioned in Sec. I, nuclear modification of F3 is not 
measured in the neutrino scattering. In hadron-hadron 
reactions, there are also issues of if-factors and final state 
interactions. In a future version of our analysis, we will 
consider to include other data. 

At this stage, the available nuclear data are mainly 
on the F2 structure functions in the electron and muon 
scattering. Experimental results are shown by the ratio 



X F 2 



M 2 ) 



F 2 A 



(x,Q 2 ) 



(x,Q 2 ) 



(3.1) 



Information about the used experimental data is given 
in Table I, where nuclear species, references, and data 
numbers are listed. The experimental data are taken 
from the publications by the European Muon Collabo- 
ration (EMC) [14-16] at the European Organization for 
Nuclear Research (CERN), the E49, E87, E139, and E140 
Collaborations [17-20] at the Stanford Linear Accelera- 
tor Center (SLAC), the Bologna-CERN-Dubna-Munich- 
Saclay (BCDMS) Collaboration [21, 22] at CERN, the 
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TABLE I: Nuclear species, references, and data numbers are 
listed for the used experimental data with Q 2 > 1 GeV 2 . 



nucleus 


experiment 


reference 


# of data 


He 


SLAC-E139 


[201 


18 




NMC-95 


[231 

L J 


17 


Li 


NMC-95 


[231 

L "J 


17 


Be 


SLAC-E139 


[201 

L J 


17 


C 


EMC-88 


[141 

L J 


9 




EMC-90 


[151 

L J 


5 




SLAC-E139 


[201 

L J 


7 




NMC-95 


[231 

L "J 


17 




FNAL-E665-95 


[251 


5 


N 


BCDMS-85 


[211 

L J 


9 


Al 


SLAC-E49 


[181 


18 




SLAC-E139 


[201 

L J 


17 


Ca 


EMC-90 


[151 


5 




NMC-95 


[231 

L J 


16 




SLAC-E139 


[20] 


7 




FNAL-E665-95 


[251 


5 


Fe 


SLAC-E87 


[171 

L J 


14 




SLAC-E140 


[191 

L 


10 




SLAC-E139 


[20] 


23 




BCDMS-87 


[22] 


10 


Cu 


EMC-93 


[16] 


19 


Ag 


SLAC-E139 


[20] 


7 


Sn 


EMC-88 


[14] 


8 


Xe 


FNAL-E665-92 


[24] 


5 


Au 


SLAC-E140 


[19] 


1 




SLAC-E139 


[20] 


18 


Pb 


FNAL-E665-95 


[25] 


5 


total 






309 



New Muon Collaboration (NMC) [23] at CERN, the E665 
Collaboration [24, 25] at the Fermi National Accelera- 
tor Laboratory (FNAL). The used data are for the nu- 
clei: helium (He), lithium (Li), beryllium (Be), carbon 
(C), nitrogen (N), aluminum (Al), calcium (Ca), iron 
(Fe), copper (Cu), silver (Ag), tin (Sn), xenon (Xe), gold 
(Au), and lead (Pb). As explained in the next subsec- 
tion, Qq=1 GeV 2 is taken as the point where the parton 
distributions are defined. Because the published data in- 
clude smaller Q 2 points, it is necessary to choose the ones 
with large enough Q 2 which could be considered in the 
perturbative QCD region. The only data with Q 2 > 1 
GeV 2 are taken into account in the x 2 analysis. The total 
number of the data is 309. 

We show the x and Q 2 points of the employed data 
in Fig. 2. The SLAC data are restricted to the large x 
and small Q 2 region. Because the SLAC data are taken 
for many nuclear species, their data are very valuable 
for our analysis. However, they cannot address the issue 
of shadowing due to the lack of small x measurements. 
The BCDMS data are also taken in the large x region. 
The difference from the SLAC measurements is that the 
BCDMS data have large Q 2 values. There exist a large 
Q 2 gap between the SLAC and BCDMS data sets, which 
may enable us to investigate nuclear Q 2 evolution. On 
the other hand, the EMC, NMC, E665 data are almost in 
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FIG. 2: Kinematical range is shown by plotting x and Q 2 
points of the used data. 



the same kinematical range. Because these data include 
small x points, it is possible to investigate the shadowing 
region as well as the medium- £ modification part. How- 
ever, the data have rather small Q 2 values in a restricted 
Q 2 range at small x. It suggests that it is difficult to de- 
termine the nuclear gluon distributions from the scaling 
violation at small x. In order to obtain the smaller x or 
larger Q 2 data than those in Fig. 2, we should wait for a 
next generation project such as HERA-eA [26] or eRHIC 
[27]. 



B. x 2 analysis 

Our analysis is done in the leading order (LO) of a s . 
The structure function F A is expressed in a parton model 

as 



F 2 A (x,Q 2 ) 



^2e 2 q x[q A (x,Q 2 )+q A (x,Q 2 )}, (3.2) 



where e q is the quark charge, and q A (q A ) is the quark 
(antiquark) distribution in the nucleus A. It is note- 
worthy to mention here that the structure functions and 
the parton distributions are defined by those per nucleon 
throughout this paper. Although it is now established 
that each antiquark distribution is different in the nu- 
cleon [28], there is no data which could suggest flavor 
dependent antiquark distributions in a nucleus [29, 30]. 
Therefore, we should inevitably assume flavor symmetric 
antiquark distributions: 



u A = d A = s A = q A . 



(3.3) 



Furthermore, the flavor number is taken three. Then, 
the structure function becomes a summation of valence- 
quark and antiquark distributions: 



F 2 A (x,Q 2 ) 



-[4u A (x,Q 2 )+d A (x,Q 2 ) + 12q A (x,Q 2 )}. 

(3.4) 
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The gluon distribution g A (x,Q 2 ) is not explicitly con- 
tained in the LO structure function; however, it con- 
tributes to F A through Q 2 evolution which is described 
by coupled integro-differential equations with g A (x,Q 2 ). 
In this way, we need four types of distributions, u A , d A , 
q A , and g A , for calculating the structure function F A . 

Because the valence-quark distributions in a nucleus 
are much different from the ones in the proton, we should 
be careful in defining the weight function Wi(x, A, Z). If 
there were no nuclear modification, in other words, if a 
nucleus were described simply by a collection of protons 
and neutrons, the parton distributions in the nucleus A 
are given by 

Atf(x, Q 2 )no-mod = Zf?(x, Q 2 ) + Nf?( X , Q 2 ), 

where fi = u v ,d v , q, or g. (3.5) 

The functions ff(x,Q 2 ) and f^ix^Q 2 ) are z-type distri- 
butions in the proton and neutron, respectively. Isospin 
symmetry is usually assumed for the parton distributions 
in the proton and the neutron: 



dy — d v , d v 



g n = ? = q, g n = g p = g. 



(3.6) 



Then, the nuclear parton distributions become 

Zu v {x,Q 2 ) +Nd v {x,Q 2 ) 



U y (x, Q )no-mod 
d v (x, Q )no-mod 



Zd v (x,Q 2 )+Nu v (x,Q 2 ) 



A 



q A (x, Q 2 )no-mod = q(x, Q 2 ), 



g A [x, Q 2 )no-mod = g(x, Q 2 ). (3.7) 

As suggested by the F A / measurements, 10—30% 
modification is expected for medium and large size nu- 
clei. The modification from the expressions in Eq. (3.7) 
should be expressed by the functions Wi(x,A,Z) at Ql 
as discussed in Sec. II: 



u A (x,Ql) = w Uv (x,A,Z) 

d A (x,Q 2 ) = w dv (x,A J Z) 

q A (x, Ql) = w q (x, A, Z) q(x, Ql), 
g A (x, Ql) = w g (x, A, Z) g(x, Ql). 



Zu v (x,Q 2 )+Nd v (x,Ql) 
A 

Zd v {x,Ql) + Nu v {x,Ql) 



(3.8) 



We would like to take small Ql value in order to ac- 
commodate experimental data as many as possible. On 
the other hand, because the Dokshitzer-Gribov-Lipatov- 
Altarelli-Parisi (DGLAP) equations are used for the Q 2 
evolution, Ql should be large enough so that perturbative 
QCD can be applied. As a point which could compromise 
these conflicting requirements, Ql =1 GeV 2 is taken. 

Next, a set of parton distributions in the nucleon is 
selected. There are available CTEQ, GRV, and MRS 
distributions. However, analytical expressions are not 



given in the CTEQ paper [1], and the GRV expressions 
are provided at small Q 2 (0.26 GeV 2 ) [2]. Therefore, we 
decided to use a LO set of MRST (Martin, Roberts, Stir- 
ling, and Thorne) [31], which is conveniently defined at 
Ql=l GeV 2 . In this way, we use the central gluon ver- 
sion of MRST-LO distributions with the scale parameter 
Alo = 0.1741 GeV in our x 2 analysis. Consequently, 
obtained nuclear distributions in Sec. IV become the 
MRST distributions in the limit A -> 1. 

Because each function Wi(x, A, Z) still has four or five 
parameters, it is necessary to reduce the total number as 
many as possible in order to become an efficient analysis. 
First, there are following obvious constraints (a), (b), and 
(c) for the nuclear distributions, 
(a) Charge 

Nuclear charge has to be the atomic number Z. It can 
be expressed in a parton model by considering an infinite 
momentum frame for the nucleus. Let us consider elastic 
scattering of a real photon with its momentum \q\ — > 
from a nucleus. Because each parton is moving very fast, 
we could use a parton picture for describing the process. 
Then, the nuclear charge is given as 



/ 
/ 



2 - {u A_ U A ) _l_ {d A_ d A ) _l_ {s A_- s A ) 



dxA 

dx ±[2u A -d A }. 



(3.9) 



The second equation is obtained because there is no net 
strangeness in an ordinary nucleus although s A {x,Q 2 ) 
could be different from 5 A (x,Q 2 ). The valence distribu- 
tions are defined by u A = u A — u A and d A = d A — d A as 
usual. 

(b) Baryon number 

Baryon number of a nucleus is A, and it is expressed 
in the parton model as 



/ 
/ 



dx A 



A 



1{U A -U A ) + ±(d A -d A ) + ±{8 A -8 A ) 



dx - [u A + d A ] 



(3.10) 



(c) Momentum 

Nuclear momentum is the addition of each parton con- 
tribution: 



A = j dxAx [u A + u A + d A + d A + s A + s A + g A ] 

= J dx Ax [u A + d A + 6 q A + g A ] . (3.11) 

If the weight functions are the quadratic functional 
type, the distributions are expressed by the parameters, 

^u v i ^dni ^vi Cy i Pv, ttqi ^qi Cq, Pq, ^gi ^gi ^Qi fig. 

Here, the valence up- and down-quark parameters are 
assumed to be the same except for a Uv and aa v because 
both weight functions are expected to be similar. How- 
ever, at least one parameter should be different in or- 
der to satisfy the conditions (a) and (b) simultaneously 
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since there are data for non-isoscalar nuclei with Z ^ N. 
Among the parameters, three of them can be fixed by the 
conditions, (a), (b), and (c). We use these three condi- 
tions for determining a Uv1 a^, and a g . 

There are still irrelevant parameters which could be 
removed from the parametrization. First, f3q and f3 g de- 
scribe the functional shape of the antiquark and gluon 
distributions at large x. However, they are irrelevant in 
the sense that both distributions are extremely small at 
large x, for example x ~ 0.8. They are not expected to 
contribute to significantly in this x region. Further- 
more, the antiquark and gluon distributions themselves 
are not determined at such large x in the nucleon. In this 
way, we decided to fix the parameters at f3q = (3 g = 1. 
We checked the sensitivity of our x 2 analysis results to 
this choice. Even if f3q = f3 g = 0.5 or 2 is taken, we found 
that the obtained x 2 changed very little. This fact indi- 
cates that x 2 is almost independent of these parameters. 
There is another irrelevant parameter. The gluon param- 
eters b g and c g determine the functional form at medium 
and large x. However, the gluon distribution does not 
contribute to F^ directly, so that the detailed x depen- 
dence cannot be determined in the analysis. Therefore, 
2c a with the following consideration. 

As far as b g 

is taken to be larger than 
the x 2 to a considerable extent. However 
smaller than — c g , it could contradict with the condition 
Wi(x — > 1, A, Z) — > +oo depending on the value of a g . In 
this way, there are seven free parameters 



b g is fixed at b g — —^ g 
The x 2 analysis tends to favor negative c ( 



c^, the choice does not affect 
if b g is taken 



bvi Cvi ftvi Q>qi bqi Cq-> C< 



■9i 



(3.12) 



in the quadratic fit. There are additional parameters in 
the cubic fit. However, as far as the gluon distribution 
is concerned, we use the same quadratic form even in 
the cubic type analysis. The structure function F^ is 
rather insensitive to the gluon distribution, especially in 
the LO analysis. It does not make much sense to in- 
troduce an additional parameter for the gluon in the x 2 
analysis without the data which could restrict the gluon 
distribution. In this way, the actual parameters are 



i Pvi ^qt bqi Cg, dq 



(3.13) 



in the cubic fit. 

In the theoretical calculations, the nuclei are assumed 
as 4 He, 7 Li, 9 Be, 12 C, 14 N, 27 A1, 40 Ca, 56 Fe, 63 Cu, 107 Ag, 
118 Sn, 131 Xe, 197 Au, and 208 Pb. The initial nuclear dis- 
tributions are provided at Q 2 = 1 GeV 2 with the param- 
eters in Eq. (3.12) or (3.13). They are evolved to the 
experimental Q 2 points by the DGLAP evolution equa- 
tions. Then, obtained structure-function ratios Rp are 
compared with the experimental values for calculating 



x 2 = E 

3 



( r)A,data 



r)A^theo\2 



(3.14) 



where the experimental error is given by systematic and 

Al- 



though the deuteron structure function is sometimes as- 
sumed the same as the one for the nucleon [6, 7], the 
deuteron modification is also taken into account in our 
analysis simply by setting A = 2. With these prepara- 
tions together with the CERN subroutine Minuit [32], 
the optimum parameter set is obtained by minimizing 

x 2 . 



IV. RESULTS 

Our analysis results are explained in Sec. IV A in com- 
parison with the used experimental data. Then, obtained 
optimum nuclear parton distributions are discussed in 
Sec. IV B. 



A. Comparison with data 

Analysis results are shown for both quadratic and cu- 
bic types. A minimal functional form is the quadratic 
type according to the discussion in Sec. II B. It is mini- 
mal in the sense that the major features of the measured 
ratios Rp 2 could be described in the whole x region in- 
cluding increase at large x ~ 0.9, depletion at medium 
at x ~ 0.6, antishadowing at x ~ 0.15, and shadowing at 
small x < 0.05 with a minimum number of parameters. 
However, there are some restrictions on the distribution 
shape, most significantly on the valence-quark distribu- 
tions. In the cubic type analysis, the shape becomes more 
flexible due to the additional parameters. Because of the 
new freedoms, obtained Xmin sn °uld be smaller for the 
cubic type, however, by sacrificing the computation time. 

Fitting results are compared with experimental data 
for all the used nuclear data in Figs. 3—16. Obtained 
optimum parton distributions are used for calculating 
the curves at Q 2 =5 GeV 2 in these figure. The dashed 
and solid curves indicate the ratios in the quadratic and 
cubic analyses, respectively. The experimental data are 




FIG. 3: Fitting results are compared with the helium data. 
The dashed and solid curves are for the quadratic and cubic 
types, respectively, at Q 2 — 5 GeV 2 , whereas the data are 
taken at various Q 2 points. 



statistical errors as ( a f ata ) 2 = (erf 8 ) 2 + (of at ) 2 
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FIG. 4: Comparison with the lithium data. 



Q 

0.9- 



a E139 


























Be 





X 

FIG. 5: Comparison with the beryllium data. 




FIG. 6: Comparison with the carbon data. 




FIG. 7: Comparison with the nitrogen data. 



taken at various Q 2 points as shown in Fig. 2, so that 
the data cannot be compared directly with the curves; 




FIG. 8: Comparison with the aluminum data. 




FIG. 9: Comparison with the calcium data. 




0.7^ 

0.001 



X 



FIG. 10: Comparison with the iron data. 



nevertheless we can see a general tendency. 

Obtained minimum x 2 values are Xmm = 583.7 and 
546.6 in the quadratic and cubic analyses, respectively, 
for the 309 total data points. Because the x 2 per degrees 
of freedom is given by x^ n /d.o./.=1.93 (quadratic) and 
1.82 (cubic), they may not seem to be excellent fits. How- 
ever, it is partly due to scattered experimental data as 
it is obvious, for example, from Fig. 9. The data from 
different experimental groups are scattered particularly 
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r <^ 0.9- 



♦ E665 


















Pb - 



FIG. 13: Comparison with the tin data. 



FIG. 16: Comparison with the lead data. 



in the small x region, and they contribute to x 2 signifi- 
cantly. Therefore, a slightly large xLn (xL/ /d -°-J P - > X ) 
is unavoidable whatever the analysis method is. 

From the figures, we find that the experimental shad- 
owing at small x < 0.05 and anti-shadowing at x ~ 0.15 
are generally well reproduced by the analysis. There are 
slight deviations from the data at medium x for carbon 
and nitrogen. However, if we try to explain the carbon 



and nitrogen data at medium x, we overestimate the de- 
pletion for the beryllium, silver, and gold as obvious from 
Figs. 5, 12, and 15. The fit is also not excellent for 
the helium at medium x. The reason could be that the 
helium nucleus is an exceptional tightly-bound system 
which cannot be explained by the simple 1/A 1 / 3 behav- 
ior. 

There are typical differences between the quadratic 
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and cubic curves in Figs. 3—16. They are different in 
the small x region, where there are not so many experi- 
mental data. The quadratic curves are above the cubic 
ones at small x (0.001 < x < 0.01), but they are below 
in the region 0.03 < x < 0.14. Both curves agree well in 
the larger x region, where there are many experimental 
data. From these figures, we find that both analyses re- 
sults are similar except for the minor differences in the 
small x region. 



B. Obtained parton distributions 

Obtained optimum weight functions in the helium, cal- 
cium, and gold nuclei are shown in Figs. 17 and 18 for 
the quadratic and cubic analyses, respectively. Because 
valence up- and down-quark functions are the same for 
isoscalar nuclei and they are very similar in other nuclei, 
only the valence up-quark functions are shown in these 





TABLE II: Each 


X 2 contribution. 


nucleus 


# of data 


X 2 (quad.) 


X 2 (cubic) 


He 


35 


55.6 


54.5 


Li 


17 


45.6 


49.2 


Be 


17 


39.7 


38.4 


C 


43 


97.8 


88.2 


N 


9 


10.5 


10.4 


Al 


35 


38.8 


41.4 


Ca 


33 


72.3 


69.7 


Fe 


57 


115.7 


92.7 


Cu 


19 


13.7 


13.6 


Ag 


7 


12.7 


11.5 


Sn 


8 


14.8 


17.7 


Xe 


5 


3.2 


2.4 


Au 


19 


55.5 


49.2 


Pb 


5 


7.9 


7.6 


total 


309 


583.7 


546.6 



He 






Ca 


Q 2 =lGeV 2 / 




-- - Au 






sea 








valence 



He 




- Ca 


Q 2 - 1 GeV 2 // 


; - - - Au 










gluon 



Each x 2 contribution is listed in Table II. We notice 
that the x 2 values per data are especially larger for the 
carbon, calcium, and gold nuclei than the average. Be- 
cause the NMC errors are very small, slight deviations 
from the NMC data produce large x 2 values. For exam- 
ple, the calcium data at x = 0.25 have peculiar behavior 
which cannot be reproduced by a smooth x dependent 
function, yet they have small errors which contribute sig- 
nificantly to the total x 2 - 

From the cubic x 2 values in Table II in comparison 
with the quadratic ones, we find significant x 2 improve- 
ments in carbon, calcium, iron, and gold, however, by 
sacrificing the x 2 values for lithium, aluminum, and tin. 
The additional degrees of freedom make it possible to ad- 
just the distribution shapes to the data. There is a 20% 
X 2 improvement in the iron from the quadratic analysis 
to the cubic one; however, the difference is not very clear 
in Fig. 10 within the region where the data exist, except 
for the x ~ 0.1 region. 

From Figs. 3—16 and Table II, we find that the anal- 
yses with the quadratic and cubic functional types are 
both successful for reproducing the major experimental 
properties. It is obvious from the x 2 comparison that 
the cubic results are better. However, because the x 2 
improvement is not very large and both curves look simi- 
lar in Figs. 3—16, both results could be taken as possible 
nuclear distributions. 



FIG. 17: Obtained weight functions for the helium, calcium, 
and gold nuclei in the quadratic analysis. Only the valence 
up-quark functions are shown as the valence distributions. 



1.2 



He // 

_ca Q 2 =lGeV 2 / 

Au / 




sea ./v V 
___________ ^7-' valence - 





0.001 0.01 0.1 1 

X 



He 




: Ca 


Q 2 - 1 GeV 2 // 


: - - - Au 










gluon 



0.001 0.01 0.1 1 

X 

FIG. 18: Obtained weight functions for the helium, calcium, 
and gold nuclei in the cubic analysis. Only the valence up- 
quark functions are shown as the valence distributions. 
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figures. 

First, the quadratic results are explained. The valence- 
quark distributions have depletion at medium x because 
they should explain the modification of the ratios Rp 
at x ~ 0.6. Because of the assumed quadratic func- 
tional form with the baryon- number constraint, the va- 
lence distributions show antishadowing property at small 
x. It indicates 2.6% antishadowing for the calcium nu- 
cleus at x = 0.001. It is noteworthy to reiterate that this 
quadratic type could not allow a shadowing property for 
the valence-quark distributions, so that it does not agree 
with a shadowing prediction, for example, in Ref. [33] 
although it could agree with a parton-model analysis of 
Refs. [6, 13] and also with the one in Ref. [34]. Next, 
the antiquark distributions should explain the shadow- 
ing of F2 at small x, so that they also have shadowing 
property at x < 0.07. The antiquark shadowing is about 
20% for the calcium. The antiquark weight functions in- 
crease as x becomes larger. They cross the line Wq = 1 
at x ~ 0.08 and continue to increase. The gluon weight 
functions have similar property to the antiquark func- 
tions except that the shadowing is smaller (7% for Ca) 
and that the crossing point is slightly larger (x ~ 0.2). 
The similar functional form is obtained partly due to the 
momentum-conservation constraint. The gluon distribu- 
tion does not contribute significantly in the LO analysis; 
however, our analysis tends to rule out the gluon anti- 
shadowing at small x. 

We performed the cubic analysis for getting better 
agreement with the data by the additional adjustable pa- 
rameters. However, as shown in Fig. 18, the obtained 
functions are similar to those in Fig. 17 except for the 
valence functions at small x. Now, the valence distribu- 
tions have freedom to have shadowing or antishadowing 
property at small x. In fact, the results show shadowing 
for the valence distributions at very small x; however, the 
magnitude is fairly small, 0.1% at x=0.001. The bump of 
Rp 2 in the region, 0.1 < x < 0.2, is explained mainly by 
the antiquark distributions. This situation is very differ- 
ent, for example, from the picture in Ref. [33], where the 
bump is explained mostly by the valence distributions. 
Of course, the F2 data are not enough to separate the 
valence and sea distributions, so that we should wait for 
future experimental activities, especially at a neutrino 
factory [9], for precise information. 

The antiquark weight functions in the cubic analysis 
are similar to the ones in the quadratic case. We ex- 
pected a possibility of much wild behavior. For example, 
it has shadowing at small x, antishadowing at x ~ 0.15, 
depletion at medium x, and rise at large x. However, 
even if the input antiquark distributions are given by 
this functional type in our x 2 analysis, they converge to 
the functions in Fig. 18. 

As mentioned in the previous section, the gluon distri- 
butions are assumed to be the quadratic functional form 
even in the "cubic" fit. Therefore, the obtained functions 
are very similar to the ones in Fig. 17. We also tried the 
cubic type for the gluon. The additional factor d g x 3 con- 



trols the behavior at large x. However, the medium and 
large x behavior is almost irrelevant for the gluon, es- 
pecially in the LO analysis of F2. Therefore, it is very 
difficult to control the gluon parameters in a meaning- 
ful way within the x 2 analysis. For example, the % 2 fit 
could produce an unphysical negative gluon distribution 
at large x if loose bounds are given for the parameters. 
It is almost meaningless to introduce the additional free- 
dom at large x without the data which could restrict the 
gluon distributions themselves. This is the reason why 
we decided to have the quadratic gluon distributions even 
in the "cubic" analysis. 

We find in Figs. 17 and 18 that the variations from the 
calcium (A = 40) to gold (A = 197) are small. Therefore, 
the obtained parton distributions could be extrapolated 
into the distributions in a nucleus with a larger mass 
number (A > 208), which is outside the analyzed nuclei 
in this paper. 

Next, errors are discussed on the obtained weight func- 
tions. As an output of the x 2 fit by the Minuit subrou- 
tine, the optimum parameters and errors are obtained. 
The error matrix has a complicated form with nondiago- 
nal elements. It is not straightforward to perform a rig- 
orous error analysis with the complicated error matrix. 
The project is, for example, in progress as an activity of 
the Asymmetry Analysis Collaboration (AAC) [4], and 
there are also recent studies in Ref. [35]. Here, we em- 
ploy a simple method which is used for example in Ref. 
[36]. Effects of the Minuit errors on wi are calculated 




0.7^ > ~ ~ I 

0.001 0.01 0.1 1 

X 




0.8- 



0.7^ > ~ ~ ' I 

0.001 0.01 0.1 1 

X 

FIG. 19: Weight functions with errors for the calcium in the 
quadratic analysis. For an isoscalar nucleus like the calcium, 
the valence up- and down-quark functions are the same. 
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FIG. 20: Weight functions with errors for the calcium in the 
cubic analysis. 

exclusively for each parameter, and then maximum varia- 
tions are shown as errors in the function Wi . The calcium 
results are shown in Figs. 19 and 20 for the quadratic and 
cubic analyses, respectively. In an isoscalar nucleus like 
the calcium, the valence up-quark function is the same 
as the valence down-quark function. The valence-quark 
functions have some errors around the minimum point at 
x ~ 0.7; however, the small x region is well determined as 
long as the functional form is fixed. However, there are 
uncertainties in the small x behavior since both valence 
functions are different at small x in Figs. 19 and 20. 
This kind of error, originating from the assumed func- 
tional form, is not taken into account in the error bands 
of these figures. 

The antiquark functions have some errors at small x\ 
however, they are obviously shadowed at small x. The er- 
rors and the distribution shapes are very similar at small 
x in Figs. 19 and 20. Other interesting point is that both 
errors are very different in the medium x region, x > 0.1. 
Therefore, the antiquark weight function cannot be well 
determined at x > 0.1, and it depends on the assumed 
functional form. There should be also differences in the 
large x region. However, because the antiquark distribu- 
tions are very small at x > 0.4 and they do not contribute 
to F2 significantly, the large-x antiquark distributions are 
not important unless we consider a reaction which is sen- 
sitive to them. 

As shown in Figs. 19 and 20, the gluon weight func- 
tions have large errors in the whole x region. The first 
reason for the large errors is that the analyses are done in 
the leading order, and the second is that only F2 data are 



used for the \ 2 analyses. Nevertheless, it is interesting to 
find that the gluon distributions are shadowed at small x 
even if the errors are taken into account. Next, there is 
a tendency of increase as x becomes larger. Determina- 
tion of large-x gluon distributions is not possible in the 
present analyses. From the simple estimate, we showed 
the errors in the weight functions. However, these studies 
are intended to give rough ideas on the errors. In future, 
we try to investigate a more complete error analysis. 

Using the results for the weight functions, we show 
the parton distributions for the calcium nucleus at Q 2 =l 
GeV 2 in Fig. 21. The dashed and solid curves are the 
quadratic and cubic analysis results, respectively. From 
the F<2 measurements, the quark distributions are rela- 
tively well determined. For determining the gluon dis- 
tributions and the details of the quark distributions, 
we need to use other reaction data. Especially, future 
hadron-collider data should be useful. 

From these analyses, we clarified how well the nuclear 
parton distributions can be determined only by the mea- 
surements of the structure functions F^. As mentioned 
in the introduction, the nuclear parametrization is still 
premature in the sense that theoretical and experimental 
efforts are necessary for determining the accurate distri- 
butions. As far as the parametrization fit is concerned, 
much detailed analyses should be done as an extension 
of our present studies. On the other hand, the authors 
hope that experimental efforts will be made for probing 
the valence-quark distributions at small x by a neutrino 
factory and for finding the antiquark and gluon distribu- 
tions by hadron colliders such as RHIC. 

Our studies should be also important for investigating 
the parton distributions in the nucleon. As mentioned 
in Sec. I, nuclear data have been used partially for ob- 
taining the parton distributions in the nucleon. In par- 
ticular, neutrino data are important for determining the 
valence-quark distributions; however, the data are taken 
for example for the iron nucleus. We need to feed back 
our studies to readjust the distributions in the "nucleon" . 



0.5- 
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FIG. 21: Parton distributions for the calcium are shown. The 
dashed and solid curves indicate the quadratic and cubic re- 
sults, respectively. For the calcium, the valence up- and down- 
quark distributions are the same. 
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V. PRACTICAL NUCLEAR PARTON 
DISTRIBUTIONS 

The nuclear parton distributions obtained in our anal- 
yses could be used for studying other high-energy nuclear 
reactions. We propose two types of distributions which 
are obtained in the quadratic and cubic type analyses. 
Because the Xmin 18 smaller for the cubic type, we pre- 
fer it to the quadratic type. The distributions are pro- 
vided either in the analytical form at Q 2 =l GeV 2 or in 
the form of computer subroutines. Although there could 
exist distributions at x > 1 in a nucleus, such large x dis- 
tributions are not provided in our studies as explained in 
Sec. II B. The analytical expressions are given in Sec. 
VA, the subroutines are explained in Sec. VB. 



A. Analytical expressions 

The analytical expressions are given at Q 2 =l GeV 2 . 
Therefore, one needs to evolve the distributions to a spe- 
cific Q 2 point by one's own evolution code with A^q ST = 
0.1741 GeV. If it were the case where Q 2 dependence can 
be neglected, one may use the analytical distributions 
without the evolution. The nuclear distributions should 
be calculated by Eq. (3.8) with the obtained weight func- 
tions and MRST-LO (central gluon) distributions [31]. 
However, one should be careful that the antiquark distri- 
butions are slightly modified from the original form so as 
to become flavor symmetric q = Sea(MRST)/6, because 
the antiquark flavor asymmetry is not taken into account 
in our nuclear analyses. 



1 Type I: Cubic fit 

We call the cubic distributions type I distributions. 
The weight functions obtained in the cubic fit are 

a Uv (A, Z) + 0.6222a; - 2.858a; 2 + 2.557a; 3 



(1 - a;) - 8107 



w dv =l+(l 



A 1 / 3 



a dv {A, Z) + 0.6222a; - 2.858a; 2 + 2.557a; 3 



(1-x) 



0.8107 



W q = 1+11- 



A 1 / 3 



-0.3313 + 6.995x - 34.17x 2 + 62.54x 3 



w a = 1 + 1 



1 — x 

1 \ a g (A, Z) + 0.8008x - 0.4004x 2 



A 1 / 3 



1-x 



The nuclear dependent constants are listed in Table III 
for the used nuclei. They depend on the mass number A 
and atomic number Z in general. 



TABLE III: Obtained parameters a Uv , dd v1 and a g for the 
used nuclei in the cubic analysis. 



nucleus 


d Uv 


a d v 


a g 


D 


-0.002178 


-0.002178 


-0.1560 


He 


-0.002178 


-0.002178 


-0.1560 


Li 


-0.002690 


-0.001716 


-0.1560 


Be 


-0.002571 


-0.001815 


-0.1560 


C 


-0.002178 


-0.002178 


-0.1560 


N 


-0.002178 


-0.002178 


-0.1560 


Al 


-0.002306 


-0.002054 


-0.1560 


Ca 


-0.002178 


-0.002178 


-0.1560 


Fe 


-0.002427 


-0.001941 


-0.1560 


Cu 


-0.002456 


-0.001916 


-0.1560 


Ag 


-0.002610 


-0.001782 


-0.1560 


Sn 


-0.002726 


-0.001686 


-0.1560 


Xe 


-0.002814 


-0.001616 


-0.1560 


Au 


-0.002902 


-0.001549 


-0.1560 


Pb 


-0.002955 


-0.001509 


-0.1560 



As obvious from the table, there are significant nuclear 
dependence in the parameters a Uv and aa v . However, the 
dependence is so small in the parameter a g that it cannot 
be shown in the table. These parameters are the same for 
isoscalar nuclei because of the conditions in Eqs. (3.9), 
(3.10), and (3.11), and this fact is clearly shown in Ap- 
pendix A. If one would like to have analytical expressions 
for a nucleus which is not listed in Table III, there are 
following two possibilities. The first method is that one 
calculates a Uv1 a^, and a g so as to satisfy the condi- 
tions of nuclear charge, baryon number, and momentum 
in Eqs. (3.9), (3.10), and (3.11) for one's chosen nucleus. 

For those who think this calculation is tedious, we pre- 
pare an alternative method. Using the three conditions, 
we find that a Uv , a^ v , and a g can be expressed in terms of 
eight integrals, which are nuclear independent, together 
with A and Z. The details of this method are so tech- 
nical that they are discussed in Appendix A. If one still 
thinks that these calculations are too much works to do, 
or if one does not have a Q 2 evolution subroutine, one 
had better use computer codes explained in Sec. VB for 
getting numerical values of the parton distributions. 



2 Type II: Quadratic fit 

We call the quadratic distributions type II distribu- 
tions. The weight functions obtained in the quadratic fit 



(5.1) 
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are 



Wo 



Wn 



1+ v 
1+6 



1 \ a„„ (A, Z) - 0.2593 x + 0.2586 x 2 



A 1 / 3 
1 

1 
1 

IV3 



(1 - X) 2108 

a d „ (A, Z) - 0.2593 x + 0.2586 x 2 
(1 - x) 2 - 108 ' 
-0.2900 + 3.774 a; - 2.236 a; 2 

a g (A,Z)+ 0.4798 a; - 0.2399 x 2 



l-x 



(5.2) 



The nuclear dependent parameters are listed in Table IV. 
We find that the A dependent variations are very small 
in these parameters of the quadratic fit. If one needs 
expressions for other nucleus, one should evaluate a Uv1 
dd v , and a g as suggested in the type I section. 



TABLE IV: Obtained parameters a Uv , 
used nuclei in the quadratic analysis. 



and a g for the 



nucleus 


a Uv 




a g 


D 


0.03745 


0.03745 


-0.09391 


He 


0.03745 


0.03745 


-0.09391 


Li 


0.03709 


0.03776 


-0.09392 


Be 


0.03717 


0.03770 


-0.09392 


C 


0.03745 


0.03745 


-0.09391 


N 


0.03745 


0.03745 


-0.09391 


Al 


0.03736 


0.03753 


-0.09391 


Ca 


0.03745 


0.03745 


-0.09391 


Fe 


0.03727 


0.03761 


-0.09391 


Cu 


0.03725 


0.03763 


-0.09391 


Ag 


0.03714 


0.03772 


-0.09392 


Sn 


0.03706 


0.03778 


-0.09392 


Xe 


0.03700 


0.03783 


-0.09393 


Au 


0.03694 


0.03788 


-0.09394 


Pb 


0.03690 


0.03790 


-0.09394 



B. Parton distribution library 



small effect, a simple linear interpolation in log Q 2 is used 
for calculating the distributions at a given Q 2 . On the 
other hand, because the x dependence is more compli- 
cated, a cubic Spline interpolation is used for calculating 
the distributions at a given x point. Running this code, 
one obtains the distributions, xu„, xd^, xq A , and xg A , 
for a specified nucleus at a given x and Q 2 point. Even 
though the antiquark distributions are flavor symmetric 
at Q 2 =l GeV 2 , they are not symmetric at different Q 2 
in the next-to- leading order [28]. However, because such 
Q 2 evolution effects do not exist in the leading order, the 
antiquark distributions are consistently flavor symmetric 
at any Q 2 . 

Second, if one would like to have the distributions in 
other nucleus, one should use the second code. Here, the 
analytical expressions in Sec. V A are used as the initial 
distributions. At first, the constants a Uv1 dd v1 and a g are 
calculated so as to satisfy the charge, baryon- number, 
and momentum conditions for a given nucleus with A 
and Z. Then, they are evolved to a requested x and Q 2 
point by the ordinary DGLAP evolution equations in Ref. 
[11]. However, one has to be careful about the requested 
nucleus in the sense that it should not be too far away 
from the used nuclei. For example, we do not support 
the distributions in an extremely unstable nucleus with 
large neutron excess. Strictly speaking, huge nuclei with 
A > 208 are also outside our supporting range. However, 
as obvious from Figs. 17 and 18, the variations of the 
parton distributions are already very small between the 
calcium with A = 40 and gold with A = 197, so that the 
extrapolation from A = 208 to nuclear matter is expected 
to be reliable. The details of the usage are explained in 
Ref. [37]. 

In the second code, it takes time for getting the results 
because the Q 2 evolution calculations consume comput- 
ing time. It does not matter to calculate the distributions 
for a few Q 2 points. However, if one would like to use it 
frequently, one may try the following. The second code 
is prepared so that one could create a grid file for a re- 
quested nucleus. Then, one can use it in the first code, 
where the computation is much faster. 



If one needs to have nuclear parton distributions in a 
numerical form at a given x and Q 2 point, one may use 
the computer codes in Ref. [37]. Two kinds of subrou- 
tines are made. First, there is a subroutine for the used 
nuclei in this paper. For other nuclei, we prepared the 
second one. 

First, if one wishes to calculate the distributions in 
the used nuclei: D, 4 He, Li, Be, C, N, Al, Ca, Fe, Cu, 
Ag, Sn, Xe, Au, and Pb, one should use the first code. In 
addition, we prepared the distributions in the nucleon be- 
cause we modified the MRST antiquark distributions as 
flavor symmetric in our studies. The kinematical ranges 
are 10~ 9 < x < 1 and 1 GeV 2 < Q 2 < 10 5 GeV 2 . The 
variables x and Q 2 are divided into small steps. Then, 
a grid data set is prepared for the parton distributions 
in each nucleus. Because the scaling violation is a rather 



VI. SUMMARY 

We have done the global analyses of existing experi- 
mental data on nuclear F2 for obtaining optimum parton 
distributions in nuclei. Assuming a simple yet reason- 
able overall A dependence, the nuclear parton distribu- 
tions are expressed in terms of a number of parameters. 
The quadratic and cubic functional forms are assumed 
for the x dependence. The parameters have been de- 
termined by the x 2 analyses. As a result, we obtained 
reasonable fit to the measured experimental data of F2. 
The valence-quark distributions are relatively well deter- 
mined except for the fact that the small x part depends 
slightly on the assumed functional form. The antiquark 
distributions are reasonably well determined at small x; 
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however, the large x behavior is not obvious from the F 2 
data. The analyses indicated that the gluon distributions 
are shadowed at small x: however, they cannot be well 
determined by the present F 2 data, especially in the large 
x region. 

We have proposed two types of nuclear parton distri- 
butions which are obtained by the quadratic and cubic 
type analyses. They are provided either by the analytical 
expressions at Q 2 =l GeV 2 or by the computer programs 
for calculating them numerically. Our analyses should 
be important not only for understanding physics mech- 
anisms of nuclear modification but also for applications 
to heavy-ion physics. Our results could also shed light 
on an issue of present parton distributions in the nucleon 
because nuclear data have been partially used in the "nu- 
cleon" analysis. 
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APPENDIX A: NUCLEAR DEPENDENT 
PARAMETERS 

From the conditions of nuclear charge, baryon number, 
and momentum, the nuclear dependent parameters, a Uv , 
dd v1 and a g , can be expressed in terms of eight integrals 
together with A and Z. It is the advantage that these 
integrals are nuclear independent. Therefore, reading the 
numerical values of the integrals and using the equations 
in this section, one can easily calculate the values of a Uv , 
dd v1 and a g for any nuclei. The necessary integrals are 
the following: 



h 



= jdx- 



H v (x) 



h 
h 



= jdx- 

-I 



(1 


— xY v 




1 


(1 






X 



u v (x), 



=Jdx- 



H v (x) 



Is = J dx ^ x y v u v(x), h = J dx — ^ x y v d v (x), 



H v {x) 



U V (X), I 6 



jdx- 



(1 


— x)& v 




1 


(1 






X 



d v (x), 



(1 - 



d v (x), 



(1 - x)Pv 

Clq+Hq(x) 



{u v (x) +d v (x)} 



1-x 



1 — x 



g(x) 



(Al) 



We should note that the parton distributions in these 
equations are those in the nucleon. The functions Hi(x) 



are the one given in Eq. (2.8^ 
numerically given in Table V. 



The integral values are 



TABLE V: 


Values of the integrals are given. 


Integral 


Type I 


Type II 


h 


-0.0007990 


-0.1474 


h 


0.008540 


-0.05297 


h 


2.406 


3.768 


h 


1.148 


1.583 


h 


0.4777 


1.157 


h 


0.1772 


0.3629 


I? 


0.08326 


-0.007650 


h 


0.5246 


0.5246 



Using these integrals, we can express the nuclear de- 
pendent parameters as 



a Uv (A,Z) = - 
a dv {A,Z) = - 
a g (A,Z) = - 



Zh + (A - Z)I 2 
Zh + (A - Z)h ' 
ZI 2 + (A- Z)h 
ZI 4 + (A - Z)h ' 
1 



a Uv (A,Z){jI 5 + (l-£)/ 6 j 



+ a dv {A,Z){jI 6 + [ 1 



h 



(A2) 



From Table V and Eq. (A2), it is possible to calculate 
the parton distributions in any nucleus. However, we 
recommend to use our results for a nucleus which is rather 
close to the analyzed nuclei. 



REFERENCES 

[1] H. L. Lai et. al. (CTEQ Collaboration), Eur. Phys. 

J. C12, 375 (2000). 
[2] M. Gliick, E. Reya, and A. Vogt, Eur. Phys. J. C5, 

461 (1998). 

[3] A. D. Martin, R. G. Roberts, W. J. Stirling, and R. 

S. Thorne, Eur. Phys. J. C14, 133 (2000). 
[4] Y. Goto et. al. (Asymmetry Analysis Collaboration 

(AAC)), Phys. Rev. D62, 034017 (2000). The AAC 

library is available at http://spin.riken.bnl.gov/aac. 

H. Kobayashi et al, error analysis in progress. 
[5] For a summary, see D. F. Geesaman, K. Saito, and 

A. W. Thomas, Ann. Rev. Nucl. Part. Sci. 45, 337 

(1995). 

[6] S. Kumano, Phys. Rev. C48, 2016 (1993); S. Ku- 
mano and K. Umekawa, SAGA-HE-130-98 (hep- 
ph/9803359). 

[7] K. J. Eskola, V. J. Kolhinen, and P. V. Ruuskanen, 
Nucl. Phys. B535, 351 (1998). 

[8] A preliminary version is reported in M. Hirai, S. 
Kumano, and M. Miyama, invited talk at the in- 
ternational symposium on nuclear physics, Mum- 
bai, India, Dec. 18-22, 2000 (see http://www- 



16 



hs.phys.saga-u.ac.jp/talkOO.html), preprint SAGA- 

HE-171-01, to be published in the proceedings. 
[9] S. Kumano, invited talk at the KEK workshop on 

neutrino factory, KEK, Japan, Sept. 13—14, 2000, 

see http:/ /www-hs. phys.saga-u.ac.jp/talk00. html 

for the presentation. 
10] I. Sick and D. Day, Phys. Lett. B274, 16 (1992). 
11] M. Miyama and S. Kumano, Comput. Phys. Com- 

mun. 94, 185 (1996). 
12] A. H. Mueller and J. Qiu, Nucl. Phys. B268, 427 

(1986). 

13] R. Kobayashi, S. Kumano, and M. Miyama, Phys. 

Lett. B354, 465 (1995). 
14] J. Ashman et al. (European Muon Collaboration 

(EMC)), Phys. Lett. B202, 603 (1988). 
15] M. Arneodo et al. (EMC), Nucl. Phys. B333, 1 

(1990). 

16] J. Ashman et al. (EMC), Z. Phys. C57, 211 (1993). 
17] A. Bodek et al. (SLAC-E87 Collaboration), Phys. 

Rev. Lett. 50, 1431 (1983). 
18] A. Bodek et al. (SLAC-E49), Phys. Rev. Lett. 51, 

534 (1983). 

19] S. Dasu et al. (SLAC-E140), Phys. Rev. Lett. 60, 
2591 (1988). 

20] J. Gomez et al. (SLAC-E139), Phys. Rev. D49, 4348 
(1994). 

21] G. Bari et al. (BCDMS Collaboration), Phys. Lett. 

163B, 282 (1985). 
22] A. C. Benvenuti et al. (BCDMS), Phys. Lett. B189, 
483 (1987). 

23] P. Amaudruz et al. (New Muon Collaboration 
(NMC)), Nucl. Phys. B441, 3 (1995); M. Arneodo 



et al. (NMC), Nucl. Phys. B441, 12 (1995). 

[24] M. R. Adams et al. (FNAL-E665 Collaboration), 
Phys. Rev. Lett. 68, 3266 (1992); H. Schellman, per- 
sonal communication (2000). 

[25] M. R. Adams et al. (E665), Z. Phys. C67, 403 
(1995). 

[26] M. Arneodo, A. Bialas, M. W. Krasny, T. Sloan, and 
M. Strikman, hep-ph/9610423. 

[27] See http://quark.phy.bnl.gov/~raju/eRHIC.html. 

[28] S. Kumano, Phys. Rep. 303, 183 (1998); J.-C. Peng 
and G. T. Garvey, hep-ph/9912370. 

[29] S. Kumano, Phys. Lett. B342, 339 (1995). 

[30] The P906 Collaboration, proposal for Drell-Yan 
measurements at the FNAL Main Injector (1999). 

[31] A. D. Martin, R. G. Roberts, W. J. Stirling, and 
R. S. Thorne, Eur. Phys. J. C4, 463 (1998). Ana- 
lytical expressions of the MRST-LO distributions at 
Q 2 =l GeV 2 are provided by R. G. Roberts and W. 
J. Stirling through personal communications. 

[32] F. James, CERN Program Library Long Writeup 
D506 (unpublished). 

[33] L. L. Frankfurt, M. I. Strikman, and S. Liuti, Phys. 
Rev. Lett. 65, 1725 (1990). 

[34] S. A. Kulagin, Nucl. Phys. A640, 435 (1998). 

[35] J. Pumplin et al, hep-ph/0101032. 

[36] Spin Muon Collaboration (SMC), B. Adeva et al, 
Phys. Rev. D58, 112002 (1998). 

[37] Nuclear parton-distribution subroutines could be ob- 
tained at the web site: http://www-hs.phys.saga- 
u. ac.jp. 



